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ABSTRACT 


A  two-dimensional  hydrodynamic  code  has  been  used  to  calculate 
the  free-air  blast  field  about  the  muzzle  of  a  5 "/51*  Kaval  gun. 

The  calculated  blast  pressure  wave  as  a  function  of  time  is  presented 
along  with  the  velocity  field.  Calculated  overpressures  and  durations 
are  compared  with  experimental  data.  The  calculations  predict  the 
formation  and  location  of  the  shock  bottle  and  formation  of  the 
"smoke  ring". 
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INTRODUCTION 


The  phenomena  associated  with  the. free-air  blast  created  by  the 
releasing  of  hot  gases  as  a  projectile  leaves  the  nuzzle  of  a  gun  has 
been  of  interest  for  some  time.  Many  experimental  investigations 
r for  example,  references  (l)  and  (2)'  have  been  undertaken  to  determine 
the  free-air  blast  parameters  such  as  overpressure  or  impulse.  These 
investigations  liave  led  to  several  empirical  formulas  for  the  free-air 
overpressure  - field  [  references-"^ 3)  and  (b)_  about  the  muzzle  of  t-  gun.. 
These  investigations  have  primarily  been  done  using  pressure  transducers , 
from  which  the  free-air  blast  parameters  of  overpressure;  positive 
duration,  impulse  and  arrival  time  of  the  shock  front  can  be  determined 
at  various  positions  about  the  gun  muzzle.  These  investigations  have 
shown  that  the  free-air  blast  field  is  two-dimensional,  that  is, 
symmetical  about  the  barrel  axis. 

Several  hydrodynamic  codes  have  in  recent  years  been  developed  for 
the.  solution  of  one  and  two  dimensional  hydrodynamic  probl»ns.  These 
codes  numerically  integrate  the- hydrodynamic  equations  of  motion,  which 
constitute  a  set  of  nonlinear  partial  differential  equations.  These 
codes  have  been  used  to  investigate  the  hydrodynamics  of  such  problems 
as  spherical  charges,  shaped  charges  and  supersonic  flows.  The  objective 
of  this  report  is  to  show  that  a  two  dimensional  hydrodynamic  code  may 
be  used  to  calculate  the  free-air  blast  field  about  the  muzzle  of  a  gun. 
Calculations  on  a  5"/5b  Naval  gun  have  been  carried  out  and  are  compared 
with  experimental  data. 

The  hydrodynamic  code  used  for  the  calculations  presented  in  this 
report  was  the  one-material  SHELL  code.  The  SHELL  family  of  hydi odyr.amic 
codes  uses  the  Eulenan  fora  of  the  equations  of  motion  which  considers 
the  mass  in  the  grid  system  to  be  a  continuum  rather  than  discrete 
particles.  Numerous  blast-wave  calculations  have  been  performed  using 
the  SHELL  codes  such  as  those  pre-ented  in  references  (5)  through  (8). 

The  SHELL  code  used  in  the  calculations  presented  here  was  obtained 
f.om  the  Air  Force  Weapons  Laboratory,  Kirtland  A.FB,  New  Mexico. 


EQUATIONS  OF  MOTION 


The  hydrodynamic  code  used,  SHELLTC,  is  a  one  material,  two- 
dimensior.al,  pure  Eulerian  code.  The  free-air  blast  about  the  muzzle 
of  a  gun  is  cylindrically  symmetric  about  the  barrel  axis  and  thus  two- 
dimensional.  The  hydrodynamic  equations  of  motion  for  axisymmetric , 
inviscid  compressible  flow  are  given  by: 
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Conservation  of  Mass 


Qlpy 

Ptnnsi 


Equation  of  State 

p  =  p  (p,I)  (5) 

where , 

p  =  density  (gm/cm3). 
u  =  radial  velocity  component  (cm/sec), 
v  =  axial  velocity  component  (cra/sec). 
r  =  radial  coordinate  (cm). 
z  =  axial  coordinate  (cm), 
p  =  pressure  (dynes /cm2). 

I  =  specific  internal  energy  (ergs/gm). 

t  =  time  (sec).  ; 

The  conservation  of  mass  is  automatically  satisfied  by  SHELLTC 
by  neither  allowing  the  creation  nor  destruction  of  mass.  Mass  which 
leaves  one  cell  enters  another  and  is  accordingly  added  to  th  e  receiver 
cell  and  subtracted  from  the  donor  cell. 
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The  axisymnetnc  grid  used  by  SHELLTC  is  shown  in  Figure  1. 


DR  and  DZ  are  the  radial  and  axial  dimensions  of  each  cell,  respectively. 
The  cross-sectional  area  of  each  cell  in  the  i"*1  column  is  given  by, 

TAU(I)  =  2-n-RC(i)-DR  (6) 

and  its  volume  by, 

VOL(I)  =  2-n-RC(i)-DR-DZ.  (7) 

The  conservation  equations  of  momentum  and  energy  are  handled  in 
two  phases  by  SHEDLTC.  The  first  phase  considers  the  fluid  at  rest  and 
determines  only  the  contributions  of  the  pressure  terns  to  the  time  deriva¬ 
tive.  The  second  phase  of  calculations  accounts  for  mass  transport.  The 
finite  difference  equations  used  in  SHELLTC  are  discussed  in  Appendix  A 
along  with  the  time  increment  calculations  and  the  stability  of  SHELLTC. 
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INITIAL  AMD  BOUNDARY  CONDITIONS 


A  starting  code,  CLAMTC,  is  used  to  set-up  the  grid  used  by  SHELLTC 
and  to  initialize  the  fluid  properties  in.-each  cell.  The  cell  size 
chosen  to  simulate  the  5"/5b  gun  blast -was  DZ  ='5.0  inches  (12;7  cm) 
and  DR  =  2.5  inches  (6.35  cm).  This  allowed  the  radius  of  the  gun  barrel 
to  be  one  cell  and  the  length  of  the  gun  tube'  to  be  5b  cells.  All  fluid 
properties  in  each  cell  were  initially-  set  to  ambient  conditions  by 
CLAMTC  including  the  interior  cells  of  the  gun  tube. 

All  boundary  conditions  and  the  interior  ballistics  at  ejection 
were  set  at  cycle  1  in  SHELLTC.  The  left  boundary  of  the  grid  (see 
-igure  A-l,  Appendix  A)  corresponding  to  the  centerline  of  the  gun 
barrel  was  set  as  reflective.*  All  other  boundaries  of  the  grid  were 
set  as  transmissive**  to  allow  flow  out  of  the  grid  through  these 
boundaries.  The  boundary  of  cells  corresponding  to  I  =  1  from  1  «_ J  £ 

5b  was  set  reflective  to  simulate  the  gun  barrel. 

The  initial  fluid  properties  of  the  gas  in  the  gun  barrel  at  time 
of  ejection  of  the  projectile  were  obtained  from  interior  ballistics  data 
from  the  5"/5b  gun  and  some  theoretical  data.  The  muzzle  pressure  at 
ejection  was  assumed  to  be  8000  psi  (5-52  x  10n  dynes/cm2).  This  value 
corresponds  with  experimental  measurements  ’reference  (9)’  of  muzzle 
pressure  at  ejection  for  the  5"/5b,  MARK  111  project  .le.  The  variation 
of  pressure  along  the  gun  tube  was  obtained  from  reference  (12).  The 
formula  for  p(x)  is: 

P(x>  =  1.0  +  Wo  fl.O  .  3A2x2  ‘ "  (8) 

e<7r  -  n  ■  6wP  \  v2  * 

>.w0  j  t 


*  A  reflective  boundary  is  one  along  which  the  normal  component  of  the 
flow  velocity  is  zero,  the  pressure  gradient  across  the  boundary  is  zero 
and  there  is  no  flux  of  the  conserved  quantities  across  the  boundary. 

**  A  transmissive  boundary  is  a  fictitious  boundary  across  which  flux 
»f  the  conserved  quantities  is  allowed  without  influence  from  points 
beyond  the  boundary. 
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p  =  pressure  at  any  point  along  bore,  lb/in2 ;■ 

T0  =  average  temperature  of  gas  at  projectile  ejection,  °R, 

R  =  gas  constant;  ft-lb/lb/°R', 

Wc  =  weight  of  propellant,  lb. 

Wp  =  weight  of  projectile,  lb, 
g  =  acceleration  of  gravity,  ft/sec2, 
i)  -  covolume,  in3/lb, 

A  =  bore  area,  in2, 

V^.  =  total  volume  of  bore  and  chamber,  in3 
x  =  distance  from  breech,  in. 

Substitution  of  the  above  quantities  for  the  5"/54  gun  into  equation  (8) 
gives  a  muzzle  pressure  of  about  9000  psi.  The  specific  impetus  at 
ejection,  RT0,  was  therefore  adjusted  to  yield  a  muzzle  pressure  correspond 
ing  to  experimental  data.  The  pressure  at  any  position  along  the  bore  is 
then  given  by 

p(x)  =  9516  -  .0208x2.  (9) 

The  gas  velocity  was  assumed  to  vary  linearly  from  zero  at  the  breech 
to  3000  ft/sec  (9.11*1*  x  106  cm/sec)  at  the  muzzle  at  time  of  ejection. 

The  gas  temperature  was  assumed  uniform  along  the  length  of  the  barrel  and 
equal  to  15260K,  which  agrees  with  experimental  data  for  the  5"/'51*  gun. 

As  noted  earlier,  a  one  material  version  of  SHELL  was  used.  The  one 
material  restraint  does  not  allow  for  the  distinction  between  propellant 
gases  and  ambient  air.  Thus  the  calculations  are  restricted  to  a  single 
equation  of  state.  The  ideal  gas  equation  of  state  with  the  ratio  of 
specific  heats  of  1.4  was  chosen  for  the  first  trial  calculations.  With 
these  assumptions,  the  mass,  density  and  specific  internal  energy  of 
each  of  the  cells  comprising  the  gas  inside  the  gun  barrel  are  calculated 
by  SHELLTC. 


After  the  initial  conditions  at  projectile^ ejection  have-been  set, 
the  flow  both  in  the  barrel  and  external,  to  the  barrel  is  determined  only 
by  the  governing  hydrodynamic  equations  and  boundary  conditions.  Thus 
both  the  flow  external  and  internal  of  the  barrel  is  calculated. 


RESULTS  AND  COMPARISON  WITH  EXPERIMENTAL  DATA 


The  numerical  calculations  tovsinulate  the  gun 'blast  from  a  5"/5*t 
gun  were  carried  out  using  the  digital  CDC  6700  computer  located  at  NHL. 
Calculations  were  carried  out  to  6.63  milliseconds  after  projectile 
ejection,  requiring  approximately  1.5  hours  of  computation 'time  on  the 
CD C  6700.  The  fluid  properties  of  each  cell  were  stored  on  magnetic 
tape  every  3r<*  cycle  (approximately  every  0.075  milliseconds  simulated 
time)  for  later  analysis  and  plotting.  The  results  of  the  calculations 
will  be  presented  in  this  section  and  compared  with  experimental  5"/5^ 
gun  blast  data. 

The  propagation  of  the  calculated  gun  blast  shock  at  k6.9°  and  90° 
from  the  line  of  fire  is  shown  in  Appendix  B,  Figures  B-l  and  B-2, 
respectively.  The  shock  fronts  are  smeared  over  3  to  U  cells  in  both 
figures  due  to  the  artificial  dissipation  in  SHELLTC.  Note  than  in 
Figure  B-l,  the  decaying  pressure  behind  the  shock  deviates  from  nearly 
exponential  at  1.3  milliseconds  (cycle  62),  giving  the  appearance  of  an 
inward  facing  second  shock.  This  is  the  initial  appearance  of  the  normal 
shock  at  the  leading  edge  of  the  shock  bottle.  At  about  5-0  milliseconds 
(cycle  210)  the  normal  shock  has  become  stationary  and  the  pressure 
signature  behind  the  shock  front  is  approaching  its  expected  shape.  The 
normal  shock  of  the  shock  bottle  appears  stationary  at  6  feet  from  the 
muzzle  or  about  lit.lt  calibers.  Experimental  measurements  [reference  (12)) 
indicate  that  the  shock  bottle  becomes  stationary  at  about  15  calibers 
from  the  muzzle.  The  oblique  shock  bounding  the  shock  bottle  near  the 
muzzle  can  be  seen  in  Figure  B-2  and  appears  about  one  foot  from  the 
muzzle.  Thus  the  hydrodynamic  code  appears  to  predict  both  the  formation 
and  location  of  the  normal  shock  at  the  front  edge  of  the  shock  bottle  and 
the  oblique  shock  bounding  the  shock  bottle. 

To  compare  the  theoretical  calculations  with  the  existing  experimental 
data,  particular  locations  about  the  nuzzle  were  chosen  and  the  calculated 
pressure  plotted  as  a  function  of  time.  Figure  B-3  shows  che  resulting 
calculated  pressure-time  history  at  5  feet  from  the  nuzzle  and  90°  from 
the  line  of  fire.  Figure  B— U  shows  the  experimentally  measured  over¬ 
pressure  f reference  (13)3  at  the  same  location.  The  calculated  peak  over¬ 
pressure  is  about  13. k  psi  and  the  positive  duration  is  1.32  milliseconds. 
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In  measuring  the  duration,  the  rise  time  of  tnevshock,  front  was  not 
included.  The  experimental  peak  overpressure  is  about  12.8  psi  and 
the  duration  is-'l.  30 -milliseconds.  Calculated' and  experimental  pressure 
traces  at  5  feet  from  the  muzzle  and  135°  from  the  line  of  fire  are 
plotted  as  Figure  B-5  and  B-6,  respectively.  The  calculat editeak  over¬ 
pressure  is  about  3.3  psi  as  compared  to  3.1  psi  experime...alij  ar.d  the 
calculated-duration  is  0.93  milliseconds  as  compared 'to  1.07  milliseconds 
measured  -experimentally. 

Figures  3-7  through  B-13  show  the  velocity  vectors  predicted  by 
SHELLTC  at  1.25,  2.0,  2.5,,  3.0,  3.5  and  U.5  milliseconds  after  projectile 
ejection,  respectively.  These- CALC0M?  plots  were -.generated  by  SHPL0T, 
a  plotting  code-develcp<;d  for  use-vith  SHELLTC.  Both  the  normal  and 
oblique  shocks  of  the  shock  bottle.-are  clearly  distinguishable  after  2.5 
milliseconds  (Figure  B-9).  The  absence  of  velocity  vectors  m  some  regions 
behind  the  gun  is  because  velocities  less  than  328  feet  per  second  t'l.O  x 
10 4  cm/sec )  are  not  plotted.  This  vas  done  because  the  size  of  the  vector 
vas  as  small  as  the  arrow  and  the  vector  magnitudes  could  not  be  dist.nguished 
in  these  regions. 

Qualitatively,  three  interesting  phenomenon  may  be  observed  in 
Figures  B-7  through  B-13.  First,  the  boundaries  of  the  shock  bottle  are 
easily  discernible  in  these  figures ,  shoving  both  the  normal  and  oblique 
shocks  described  earlier.  Second,  the  shock  front  is  shovn,  and  also  its 
propagation  vith  time.  Third,  a  circulating  flov  "smoke  ring"  is  shewn 
which  is  observed  in  many  gun  firings. 

Quantitatively,  few  if  any  velocity  measurements  of  the  gas  outside 
the  gun  barrel  have  been  made.  Thus,  no  experimental  data  exists  to 
compare  vith  the  calculated  velocity  field.  The  calculations  do,  however, 
predict  the  observed  phenomena  discussed  above. 

Figures  B-lt  and  B-15  show  the  calculated  and  experimentally 
measured  shock  front  positions,  respectively,  as  a  function  of  time.  The 
calculated  shock  is  seen  to  be  somewhat  slower  than  the  actual  shock, 
particularly  aft  of  the  muzzle.  The  general  difference  in  shock  velocity 
may  be  attributed  to  the  simple  equation  of  state  used  in  SHELLTC  and 
the  limitation  to  one  material.  The  larger  difference  aft  of  the  muzzle 
may  be  due  to  the  absence  of  che  effects  of  the  projectile  on  the  shock. 

T.  D.  Taylor  reference  (lit)’  has  numerically  solved  the  inviscid  flow 
equations  and  calculated  the  muzzle  blast  for  a  It. 2  inch  mortar  with  and 
without  a  shell  m  the  flow.  His  results  indicate  a  stronger  blast  field 
aft  of  the  muzzle  with  a  shell  than  without  and  thus  a  larger  shock 
velocity  behind  the  muzzle.  However,  no  comparisons  are  made  with 
experimental  data  in  reference  (lU). 
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The  stability  of  SHELLTC-  is  diseased  in. Appendix  A.  The  energy 
dissipating  parameters  SI,  S2  and  S3  (see  Appendix  A)-  were'-all  set  to 
zero  for  the- calculations 'presented  in  this  report.  Thus  the  only 
dissipation  effects  were  those  inherent  in  SHELLTC.-  Close  examination 
cf  Figure  B-10  shows  that  one  velocity  vector  near  the  axis  of  symmetry 
is  negative.  This  is  a. result  o*f  an  instability  which  occurred  at  about 
cycle  125  and  can  be  seen  to,  grow  with  time  from  Figures  B-li  through 
B-13.  The  instability  appears  to  grow  as  a, disturbance  in  the  flow. 

It  did  r.ot  effect  the  shock  front, 'at  least  out  to  cycle  260,  which- was 
of  primary  interest  in  these  calculations.  No  attempt  was  made  to 
resolve  t  -e  instability  because  of  the  lengthy  running  time  on  the 
computer.  However,  variation  of  the  energy  dissipating  parameters  to 
increase  the  artificial  dissipation>may  resolve  the  instability. 


CONCLUSIONS  AND  PLANS  FOR  FUTURE-WORK 


The  two  dimensional  hydrodynamic  code,  SHELLTC,  has  been  shown 
suitable  for  the  calculation  of  a  free-air  muzzle  blast.  Reasonably 
good  agreement  was  obtained  between  calculated  and  experimental  values 
of  peak  overpressure  and  duration  at  several  locations  about  the  muzzle. 

The  initial  conditions  in  the  barrel  at  projectile  ejection  were  obtained 
from  experimental  measurements  and  interior  ballistics  theory.  Because 
the  full  length  of  the  gun  barrel  was  simulated  by  the  grid;  no  assumptions 
were  necessary  as  to  the  variation  with  time  of  the  fluid  properties  in  the 
gun  barrel. 

The  actual  cost  of  computer  time  tc  simulate  the  5"/5**  gun  blast 
presented  in  this  report  was  about  $1200  on  the  CDC  6700  computer.  This 
cost  is  not  unreasonable  compared  to  firing  one  round  from  the  5",'.- 
gui.  To  obtain  similar  information  experimentally  would  require  multiple 
channels  of  instrumentation  and  still  would  not  produce  but  a  small  portion 
of  the  data  obtained  from  the  calculations.  The  finite  difference  cal¬ 
culations  easily  allow  the  variation  of  parameters  such  as  muzzle  pressure, 
gas  velocity,  thermal  energy  and  barrel  length,  to  determine  their  effects 
on  the  gun  blast  field.  Also,  the  computer  provides  data  on  the  velocity 
field,  which  cannot  be  obtained  using  conventional  instrumentation. 

Additional  calculations  using  a  multi-material  version  of  SHELL  are 
planned.  These  calculations  will  include  more  precise  equations  of  state 
for  each  material  and  a! low  chemical  reactions  to  occur  between  materials. 
These  calculations  should  predict  the  frec-air  blast  field  more  precisely 
and  give  further  insight  of  the  shock  bottle  and  secondary  blast,  due  to 
muzzle  flash. 
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There  are  many  applications  for  the  predicted  gun  blast  field. 

In  addition -to  gaining  insight  to  the  blast  phenomena,  propellant  gas 
influences  on  the  projectile; outside  the  muzzle,  may  now  be  calculated. 
The  predicted  blast  field  presented  in  this  report  has  already  been 
used  to  predict  the  effects  of  the  propellant  gas  on, the  motion  of 
fragments  of  a  rotating  band  which  has  been  torn  from  a  projectile  at 
ejection  from  the  muzzle. 
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APPENDIX  As 


FINITE  DIFFERENCE  EQUATIONS  OF  SHELLTC 


The  conservation  equations  are  handled  in  two  phases  by  SHELLTC. 

The  first  phase  considers  the  fluid  at  rest  and  determines  only  the 
contributions  of  the  pressure  terms  to  the  time  derivative.  Dropping 
the  transport  terms  in  .the  momentum  and  energy  equations  see  equations 
(2),  (3)  and  (!*)  ‘  ,  they  become,  in  finite  difference  form,  the  following: 

Momentum 


U  =  U(n)  +  1  -  -  PR<n>  ^  Ka-1) 

K  K  p(n)  '  dr 


v  „  v(»)  +  1  At  ( A— S) 

K  K  7^  AZ 


Energy 


I  =  I  +  Pk  '  -  vrt”)  +  VB-VT 

K  K  (n)  At 


+  2  Ul/n)  -  URtn)  +  ut-ua  •  m 

ri  +  rj_^  dr  2 


where 


PL  = 

PK 

+ 

PKL 

2 

PB  = 

PK 

+ 

PKB 

2 

VB  = 

VK 

+ 

VKB 

2 


PR  = 

PK 

+ 

PKR 

2 

FT  = 

PK 

+ 

PKT 

2 

VT  = 

+ 

VKT 

2 


(A-3) 


A-l 


UL  =  Vri  +  ri-l)  +  Um(  ri -1  +  ri-2> 

j. 

UR  =  *  ri_i)  +  UKL(ri+1  +  ri) 


K  is  the  index  of  the  cell  center  and  i  is  the  index  of  the  right-hand 
boundary  of  the  Kt*1  cell.  The  adjacent  cells  to  the  Xth  cell  are 
identified  in  Figure  A-l. 

Two  successive  passes  are  made  through  the  first  phase  of  calculations 
jf  each  cell.  New  velocities  are  calculated  from  the  momentum  equations 
for  a  full  time  step.  The  internal  energy  is  calculated  for  a  half  time 
step  using  the  old  velocities.  The- internal  energy  is  then  calculated 
for  the  remaining  half  time  step  using  the  new  velocities.  After  each 
pass  through  phase  one,  the  internal  energy  of  the  cell  is  checked  to 
insure  that  it  has  not  become  negative.  Should  such  an  event  occur,  a 
negative  time  step  is  used  to  return  the  state  of  the  fluid  to  its 
original  state  at  the  beginning  of  the  cycle  and  a  smaller  time  step  is 
calculated. 

The  second  phase  of  calculations  accounts  for  mass  transport.  The 
equation  for  conservation  of  mass  [equation  ,(l)j  becomes,  in  finite 
difference  form. 


&S.  =  '  ri-l  Pi-1  ul-l  _  ri  Pi  "i  +  P.1-1  v.1-l  P.1v,1  ' 

At  • -  -  - -  * 


(A-l*) 


A  r 


ri-h 


At 


The  conservation  of  mass  equation  can  be  written  in  another  form  as 


4n~x)  =  4n)  * 


2  2  —  — 

HAt  i  (r,  -  r.  ,)  (V  p  -  v  P  ) 

u  i  i-1  B  B  T  T 


+  2Az  (ri_1  u  LoL  -  rt  U  RpR>] 


(A-5) 


A-2 


-  i  ,  '«■  '■  yt' 

■r ,  aaiaua  ajfeaaata 


R 


FIGURE  A-l  SHELLTC  Mesh 


where  vp  'is  the  mass  flux  across  a  face  of  cell  k  and  the  subscripts 
B,  T,  L  and  R  designate  bottom,  top,,  left  and  right  faces ,  respectively. 
The  densities  are  those  of  the  donor  cells  and  the  velocities  are 'given 
by: 


7  =  1/2  (vT  * 

1  +  ^VR  “  vk'|| 


-  _  1/2  (uR  +  Uj^) 


1  +  (u  -  u  )M 
R  K  Az 


(A -6) 


This  velocity  weighting  scheme  is  used  to  ensure  stability  in  regions 
behind  a  shock  front. 

Calculation  of  the  mesh  cells  by  SHEU.TC  proceeds  from  bottom  to 
+  ":p  and  from  left  to  right.  Thus  mass  transport  for  each  cell  is  calculated 
for  the  top  and  right  faces  only.  The  new  mas3  in  the  Kth  is  given  by 


(n+l)  (n) 

Mk  =  Mjj  +  AMg  +  AMg  -  AIAp  -  AMg  ( A— T ) 


where  the  AM's  are  the  masses  transported  across  the  respective  faces  of 
the  Xtn  cell. 

New  velocities  are  computed  using  the  conservation  equations  for  axial 
and  radial  momentum  as  follows: 


v(n-H)  =  MK<n)  vK  +  AMgVg  +  AMjVg  -  AMj  vT  -  AMRvR 
M(n+1) 


(A-8) 


A-l» 


The  new  specific  internal  energy  is  computed  by  adding  the  total 
energy  carried  by  the  transported  masses  to  the  kth  cell  minus  the  new 
kinetic  energy  and  then  dividing, by  the  new  mass. 


(n+1) 


MK(n)  ilK  +  +  vjp]  +  ilig  {lB  +%(>4  +  vb). 


+  At^  [  IL  +  H  (uB  +  vB)]  -  Af^,  { IT  +  h  (u^  -  v^)j 


^R  iXR  +  *  (uR  +  VR)'}// 


(n+1) 


.  h  fu(n+1)  l2  -  ^ ■  •M+1'>  f  \ 


(n+1) 


(A-10) 


The  time  step,  it,  is  determined  by  consideration  of  stability  to 
satisfy  two  conditions.  First,  the  Courant  condition  which  prohibits  the 
transmission  of  a  signal  across  more  than  one  cell  in  one  time  step. 
Second,  that  |u/Ar|max  and  |v/Az [„ax  be  less  than  1/At,  thus  prohibiting 
the  transport  of  mass  across  more  than  one  cell  in  one  tine  step. 


STABILITY 


Von  Neumann  and  Richtmeyer,  reference  (10),  have  shown  that  the 
hydrodynamic  equations  can  be  solved  numerically  if  an  artificial 
dissipation  tern  is  introduced.  This  in  effect  smears  the  shock  front 
making  the  fluid  properties  such  as  pressure,  density,  entropy  and 
internal  energy  continuous  across  +,he  shock.  Thus  the  Rankine-Hugoniot 
conditions  are  satisfied  across  the  shock  and  the  numerical  calculations 
proce  -d  as  if  no  shock  were  present. 

The  first  phase  of  SHELLTC  calculations  are  unstable,  because  they 
contain  no  dissipative  mechanism.  However,  the  second  phase,  which  deals 
with  the  mass  movement,  has  been  shown  by  Harlow,  reference  (11),  to  give 
stability  to  the  calculations  by  adding  effective  viscosity  and  effective 
heat  conduction. 


In  addition  to  the  inherent  dissipative  effect  in  SHELLTC,  Havrocki, 
i eference  (5),  has  added  an  artificial  viscosity  scheme  to  SHELLTC.  By 
appropriate  selections  of  the  three  parameters  SI,  S2,  and  S3, 'the  energy 
dissipating  effect  may  be  varied  to  .enhance  or  diminish ”ihe  inherent 
effective  dissipation  in  SHELLTC. 

Leigh,  reference  (8),  has  conducted  a  parametric  study,- of  a  one¬ 
dimensional  planar  blast  wave  to  determine  the  effect  of  variations  in 
these  parameters.  Kis  results  are  presented  in  reference  (8),.  A  listing 
of  the  basic  SHELLTC  code  is  also  given  in  reference  (8). 
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